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3 . 

, An approach is developed for the determination of the current flowing through a nanosize silicon- 

on-insulator (SOI) metal-oxide-semiconductor field-effect transistors (MOSFET). The quantum me- 
l/") ' chanical features of the electron transport are extracted from the numerical solution of the quantum 

Liouville equation in the Wigner function representation. Accounting for electron scattering due to 
ionized impurities, acoustic phonons and surface roughness at the Si/SiC>2 interface, device charac- 
teristics are obtained as a function of a channel length. From the Wigner function distributions, the 
coexistence of the diffusive and the ballistic transport naturally emerges. It is shown that the scat- 
tering mechanisms tend to reduce the ballistic component of the transport. The ballistic component 
increases with decreasing the channel length. 
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I. INTRODUCTION 



Down scaling MOSFETs to their limiting sizes is a key challenge of the semiconductor industry. Detailed simu- 
lations that capture the transport and the quantum mechanical effects that occur in these devices can complement 
' experimental work in addressing this challenge. Furthermore, a conceptual view of the nanoscale transistor is needed 
to support the interpretation of the simulations and experimental data as well as to guide further experimental work. 
The objective of this work is to provide such a view by formulating a detailed quantum-mechanical transport model 
O , and performing extensive numerical simulations. We develop a model along these lines for the nanosize SOI, or 
thin film MOSFET, since the suppression of the short-channel effect seems to be most adequately addressed by this 
y—{ [ technology. Therefore, understanding the SOI MOSFET including its quantum-mechanical transport properties is 
>■ urgently needed. In this work, we restrict our attention to the steady-state current-voltage characteristics. 

The SOI MOSFET^ 1 is a transistor, built on a thin silicon layer, which is separated from the substrate by an 
insulating layer (Fig. 1). In the nanoscale transistor, we have to optimize device functionality and reliability. To 
achieve this goal, we need to suppress any short-channel effects as much as possible. These effects include threshold 
voltage variations versus channel length and the drain-induced barrier lowering effect, i.e. a drain voltage dependence 
■ of threshold voltage, which complicates transistor design at a circuit level. Frank et alSk recently showed, that 
increasing the channel doping concentration, -/Vch, suppresses short-channel effects. The integration of the down- 
scaled devices requires that the gate insulating layer thickness should be reduced and/or higher dielectric constant of 
the insulator are implemented. 

Si-based MOSFETs with typical sizes about 100 nm have found an application in highly integrated systems. 
Mechanisms of the electron transport in these devices differ from those in the devices of 50 nm and below. The 
'conventional' devices are described by the Boltzmann transport equation and approximation thereof. This theories 
focus on scattering-dominant transport, which typically occurs in long channel devices. On the contrary, in a structure 
£j \ with a characteristic size of the order of the mean free path, the electron transport is essentially quasi-ballistic. 

To describe the electron transport through a nanoscale transistor, a quantum-mechanical treatment is required. 
This includes the following approaches: the non-equilibrium Green's function formalism of Keldysh and Kadanoff- 
Baym&£i£i§iL2i2ii2iAiii£, the Pauli master equation methocUiii 4 * 1 ^, the full density matrix methoojiSiii and Wigner 
function method 18 i 19 i 20 . 21 i 22 i 23 i 24 . 

The main goal of the present work is the investigation of quantum transport in a nanosize MOSFET. Consequently, 
quantum effects such as quantum reflection and quantum tunnelling, occur in the electron transport. Simultaneously, 
electron scattering should be taken into account. Due to a low electron-scattering rate, the sub-50 nm devices are 
expected to offer higher speed performance than the conventional devices do. In a conventional MOSFET, the 
conductance is determined by the scattering rate of electrons in the channel, and the velocity of charge carriers is 
limited to a maximal value of about 10 7 cm/s. In MOSFETs with the channel length less than the mean free path, 
a velocity overshoot effect occurs^. The small length of the channel results into an increased on-state current and 
a reduced gate capacitance. For example, in Ref.— the fabrication of sub-0.1 /mi MOSFETs was reported. The 
resulting devices demonstrate a rather good performance in a wide range of temperatures. 
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In 2 *, a flexible two-dimensional model has been developed, that optimally combines analytical and numerical meth- 
ods and describes the key features of cylindrical SOI MOSFET devices. The theoretical tools to describe the quantum 
transport features involve the Wigner distribution function formalism 2 ^, in which a Boltzmann collision term represent- 
ing acoustics phonons, impurity and surface-roughness scattering was incorporated. Therefore this model is capable 
to probe the competition between three effects: quantum reflection, quantum tunnelling and phase de-coherence due 
to elastic scattering. The Wigner function formalism offers many advantages for quantum modelling 27 - 28 . First, it is a 
phase-space distribution, similar to classical distributions. Because of the phase-space nature of the distribution, it is 
conceptually possible to use the correspondence principle in order to determine, where quantum corrections enter the 
problem. At the boundaries, the phase-space description allows for separation of incoming and outgoing components 
of the distribution, which therefore permits the modelling an ideal contact, and hence an open system22. By coupling 
the equation for the Wigner function to the Poisson equation, we obtain a fully-self-consistent model of the SOI 
MOSFET. 

This paper is organized as follows : in Section II, a description of the system is presented in terms of a one-electron 
Hamiltonian. In Section III the quantum Liouville equation satisfied by the electron density matrix is transformed 
into a set of one-dimensional equations for partial Wigner functions. The boundary conditions for the Wigner function 
are considered in Section IV. The ballistic regime of the structure described in Section V. A one-dimensional collision 
term is derived in Section VI, in which we extend earlier work24 by also including scattering due to interface roughness. 
In Section VII, we describe a numerical model to solve the equations for the Wigner function. In Section VIII, the 
results of the numerical calculations are discussed. 



II. THE HAMILTONIAN OF THE SYSTEM 



We consider a n-channel SOI MOSFET structure. When a positive gate voltage is applied, electrons in the channel 
are confined to a narrow inversion layer near Si/Si02 interface. The current flows in the z-direction, the confinement 
is in the ^-direction, and the width of the transistor is found along the y-direction. The source and drain regions with 
the lengths Ls and Ld , respectively, are n-doped by phosphorus with concentration Nu , whereas the p-doped channel 
with the length Lea has concentration N A of boron acceptors. The lateral surface of the semiconductor brick in the 
channel region is covered by the Si02 oxide layer, while the aluminum gate overlays the oxide layer. 

In the semiconductor, the electron motion is determined by the following Hamiltonian : 
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where r = (x, y, z); V(r) = Ve(r)-|-I4(r) is the potential energy associated with the energy barrier and the electrostatic 
field; rn* x , mj y and m™ are the effective masses of the motion along the x, y and z-axes, respectively, of an electron 
of the j-th valley. 

The Schrodinger equation is solved within an effective mass approximation. It is assumed that Si/Si02 interface 
is parallel to the (100) plane. The conduction band in bulk silicon can be represented by six equivalent ellipsoids. 
When an electric field is applied in the [100] -direction, these six equivalent minima split into two sets of subbands^. 
The first set of subbands (unprimed) is two fold degenerate and represents those ellipsoids that respond with a heavy 
effective mass in the direction of the applied electric field. Because of the heavier mass, the 'unprimed' subbands have 
relatively lower bound-state energies, as compared to the 'primed' subbands and are therefore primarily occupied by 
electrons^. 

The electrostatic potential energy V e (r) satisfies Poisson's equation 

V.[ ei W c (r)] - iY-»(r)+JV D (r)- N A (r)\ , i = l,2, (2) 

where E\ and £2 are the dielectric constants of the semiconductor and oxide layers, respectively; n(r), A^D(r)and Na(t) 
are the concentrations of electrons, donors and acceptors, respectively. The barrier potential is non-zero in both oxide 
layers, and its value is denoted by a constant Vb- In our calculations, we assume that the source electrode is grounded 
whereas the potentials at the drain and gate electrodes equal Vb and Vq respectively. The system is assumed to be 
uniform in the y-direction. Consequently, the electron density is constant along the y-axis and the edge effects in that 
direction are supposed to be negligible. 

The study of the charge distribution in the cylindrical MOSFET with closed gate electrode in the state of the 
thcrmodynamical equilibrium (see Ref^i) has shown that the concentration of holes is much lower than that of 
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electrons, so that electron transport is found to provide the main contribution to the current flowing through the SOI 
MOSFET. For that reason, holes are neglected in the present transport calculations. 



III. THE QUANTUM LIOUVILLE EQUATION 

In order to study the device architecture starting from the Wigner function formalism, an equation is needed that 
describes the response of the Wigner function to changing external conditions. The time evolution of the Wigner 
function is derived from the quantum Liouville equation by applying the Wigner- Weyl transformation. Neglecting 
intcr-valley transitions in the conduction band, the one-electron density matrix can be written as 

p(r,r') = ]>>(r,r'), (3) 

3 

where Pj(r, r') is the density matrix of electrons residing in the j-th valley satisfying the quantum Liouville equation 



(4) 



In order to impose the boundary conditions for the density matrix at the electrodes, it is convenient to describe the 
quantum transport along the z-axis in a phase-space representation. For that purpose, we expand the j'-valley density 
matrix into a series with respect to a complete set of orthonormal functions ^j pa (x, y, z) and rewrite the density 
matrix pj in terms of £ — (z + z')/2 and rj = z — z' coordinates as 

1 f +oa 

Pj( r , r ') = V / dk ex P( [kr l) fjpsp's>((,k)t>jp S (x,y,z)-$* jp , s ,(x',y',z') . (5) 

psp s 

According to the symmetry of the system, these functions take the following form : 

$j P s(x, V, z) = — L= ipj s (x, z) exp (ipy) . (6) 
The functions ipj S (x, z) are chosen to satisfy the equation 

h 2 d 2 
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ipj S (x, z) + V(r, z)ip JS (r, z) = Ej a (z)il>j a (x, z) (7) 



that describes the confined motion of an electron. Here Ej S (z) are the eigenvalues of Eq. (JJJ for a given value of 
the z-coordinate that appears as a parameter. It will be shown, that Ej S (z) plays the role of an effective potential 
in the channel, and that ^>j ps (x, y, z) is the corresponding wavefunction of the motion along the a;-axis at a fixed z. 
Substituting the expansion (0 into Eq. and using Eq. Q, we arrive at the equation for fj psp > s >((,k) : 

df s r s r (C,k) hk d 1 f , , , 

gl 1 = ~^^~Q^hpsp's'{(i k) + ^ J Wj pap > a >((,k- k)fj pap i a >((,k)dk 

— oo 

+oo 

- E / Ki's' (c. k, *o Wsi (c, k') dk' , (8) 



where 



1 



W jpap , a ,((, k-kT) = — I d V (E jpa (C + Ty/2) - E jp , a ,(( - r,/2)) exp (i(fc' - k)rf) , (9) 



1 
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■Tj/2,k') 



r)/2, k') exp (i(fc' — k)rf) drj, 



(10) 
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and 



t pssi (C + v/2, fc') = ^rb 3SSl (C + 77/2) + ^p SSl (C + ry/2) (-i^ + 2 fc') , (11) 
t; s , s , ( C - „/2, fc') = 2^j-^ Sl (C - 17/2) + ^c JSSl (C - r?/2) (i^ + 2tf) 



Q2 

bj SSl {z)= I ip* s (x lZ )— ip jsi (x,z)dx , (12) 



Cj 5Sl (z)= / i/)* 8 (x,z)g- tp jSl (x,z)dx , (13) 



£ ips (z) = E js {z) + . (14) 

The first term in the right-hand side of Eq. © is derived from the kinetic-energy operator of the motion along the 
z-axis. It is exactly the same as the drift term of the Boltzmann equation. The second component plays the same role 
as the force term does in the Boltzmann equation. The last term in the right-hand side of Eq. © contains the kernel 

Mjpspi s /(C,k,k'), which mixes the functions fjp Sp > s ' with different indices s,s'. This term appears because ^pj S (x,z) 
are not eigenfunctions of the Hamiltonian 

IV. BOUNDARY CONDITIONS 

To describe the behavior of the SOI MOSFET through solving Eq. ©, we need to specify the boundary conditions 
for the functions fj pS p's'((, fc) that permit particles to enter and leave the system. The SOI MOSFET is described as 
an open system. Being part of an electrical circuit, it exchanges electrons with the circuit. For the present purposes, 
the term "open system" is used here to characterize a system that is connected to contacts (reservoirs of particles) and 
the interaction between the system and a contact necessarily involves a particle current through an interface between 
the system and contacts. 

The quantum Liouville equation (in the Wigner-Weyl representation) is of the first order with respect to the 
coordinate £ and does not contain derivatives with respect to the momentum. The characteristics of the derivative 
term are lines of constant momentum, and we need to supply one and only one boundary value at some point on 
each characteristic, because the equation is of first order with respect to the coordinate Thus, the Wigner function 
is a natural representation for an open system (SOI MOSFET). The implementation of the above described picture 
and the comparison of Eq. JSJ with the corresponding expansion of the density matrix in the equilibrium state, leads 
to the following boundary conditionsi2i22i22: 

fjmsm's'(0, fc) = Zdss'Smm' [exp P{Ej sph - E FS ) + l] 1 , fc > , 

fjmsm's'(L, fc) = 25 sa >6 mm > [exp f3(E Jspk - E FB ) + l] 1 , fc < , (15) 

where the total energy is Ej sp k = /i 2 fc 2 /2m| z + ft 2 p 2 /2mj y + Ej s (0) for an electron entering from the source electrode 
(fc > 0) and Ej S pk — h 2 k 2 /2m| z + ft 2 p 2 /2mj y + Ej a {L) for an electron entering from the drain electrode (fc < 0). In 
Eq. (|15fl (3 = 1/fceT is the inverse thermal energy, while £fs and £fd are the Fermi energy levels in the source and 
in the drain respectively. Note, that Eq. (|15f) meets the requirement of imposing only one boundary condition on the 
function fj pS p's'{(, fc) at a fixed value of fc as Eq. ((SJ) is a first-order differential equation with respect to £. 

V. BALLISTIC REGIME OF THE SOI MOSFET 

The functions fjp S p> s '(Ci k), which are introduced in Eq. (jSJ), are used in calculations of the current and the electron 
density. The expression for the electron density follows directly from the density matrix as n(r) = p(r,r). In terms 
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of the functions fj pS p's'{Ci k), the electron density can be written as follows : 

+00 

n ( r ) = ^~ / fipsp's'( z > k ) dk ^jps(x,y,z)^j pls ,(x,y,z) . (16) 

The current density can be expressed in terms of the density matri»S : 



2m.,- V <9r <9r' 



The total current passing through the cross-section of the structure at a point z, can be obtained by an integration 
over the transverse coordinates. Substituting the expansion (J5J into Eq. I|17l) and integrating over x and y, we find 

+00 +00 

J = e J2^t J dk ^™f jpsps ( Z,k ^~^™ ^ c JSS ,{z) J dk %(f jpsps >{z,k)) , (18) 

j.p.s _ oc i 3 j,p,s<s' _ QO 

where 3(/) is the imaginary part of /. The first term in the right-hand side of Eq. I|18[l is similar to the expression 
for a current of the classical theory^. The second term, which only depends on the non-diagonal functions fj psp s' , 
takes into account the effects of intermixing between different states of the transverse motion. 

Since M s ^ s s ^, s , couples functions fjpsp's'iQ, k) with different quantum numbers (jps), it can be interpreted as a collision 
operator, which describes transitions of electrons between different quantum states of the transverse motion. As shown 
in Ref^i, the third term in the right-hand side of Eq. ijSJ is significant only in the close vicinity of the p-n junctions. 
Therefore, this term is assumed to give a small contribution to the charge and current densities. Using the above 
assumption, we can treat the last term in the right-hand side of Eq. |JSj as a perturbation. Here after, we investigate 

the steady state of the system in a zeroth order approximation with respect to the operator Mj^, s ,. Neglecting the 
latter, one finds that, due to the boundary conditions (|15|) . all non-diagonal functions fj pS p>s'((, k) (p p' or s ^ s') 
are equal to zero24. Furthermore, summation over p in Eq. (JSJ gives 

+00 

^^/i.(C,fe) ~\J W js (C,k-k>)f js (t,k') dk' = (19) 

— OO 

with 

fjs(C,k) = ^J2f3Ps P s(Ck). (20) 

V 

In Eq. I|19(l the following notation is used : 

+00 

Wjs(t,k)=-^ J (^(C + ry/2)-^(C-V2))sin(fe 77 )dr 7 . (21) 

—00 

The effective potential Ej S (z) can be interpreted as the bottom of the subband (j,s) in the channel. The function 
fjs(C> k) is referred to as a partial Wigner function describing electrons, that are traveling through the channel in the 
inversion layer subband (j, s). 

VI. ELECTRON SCATTERING 

In this section, we consider the electron scattering due to acoustic phonons, impurities and surface roughness at 
the Si-SiC>2 interface. For this purpose, we introduce a Boltzmann-like single collision termSiM, which in the present 
case has the following form 

St fjsmk ^ ^ (Pjspk.j' s'p'k' fj' s'p'k' Pj spk,j' s' p' k' f j spk) • (22) 

j' s'p' k' 
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As noted above, we have neglected all transitions between quantum states with different sets of quantum numbers j 
and s. In general, the density matrix has off-diagonal elements between different subbands in the presence of inter- 
subband scattering processes. But when the broadening of inverse layer levels due to inter -subband scattering is 
sufficiently small in comparison with subband energy separations, we can neglect such off -diagonal parts and use the 
transport equation (|19() . which becomes a set of decoupled equations for distribution functions associated with each 
subband (j, s), i.e. in this case the electrons in the different subbands can be considered independent current carriers. 

The transition rates depend upon the full set of quantum numbers of each three-dimensional state. Because the 
present calculations consider only the longitudinal momentum fc, the scattering rates must be "projected" onto the 
one-dimensional model. To do so, we first assume that the distribution of electrons with respect to the motion along 
the y-dircction is the normalized Maxwellian distribution at a fixed temperature, because in the source and drain 
contacts the distribution of electrons over the quantum states of the motion along the y-direction corresponds to an 
equilibrium distribution. Consequently, due to the symmetry of the system, it is reasonable to assume that across the 
whole structure the electron distribution is given by 



fjspk(C) = fjs((,k) w 



jsp , 



(23) 



where the integration measure 



/ h 2 i3 

2mJ y 7r 



exp 



f3h 2 P 2 

2mY 



dp 



(24) 



is the normalized Maxwellian distribution function with respect to the momentum p. The integration of Eq. (|22|l over 
the angular momentum gives the collision term 

St f js (C,k) =J2(Pjs(Ck,k')f js (C,k') - P js (C,k',k)f js (C,k)) , (25) 



where 



VP 



Pjs(.Ci kyk^) — ^ ^ Pjspkjsp' k'lVjsp' ( 



(26) 



This collision term is directly incorporated into the quantum Liouville equation l|19|) as 

W js (C, k, k') = W js (C, k-k') + ^ (p ja ({, k, k') - 5 k , k , P js (C, k', k)j 



(27) 



where Wj S (z, k, k') is the modified force term in Eq. 1(19(1 . 

The scattering rates in the first Born approximation^ are evaluated according to the Fermi's golden rule 



P 



jsmkjsm' k' 



2lT ~ 2 

— {jsp'k'\H int \jspk} S (Ej B p'k' - Ejspk) , 



(28) 



where i?i n t is the Hamiltonian of the electron-phonon, the electron-impurity or electron-interface roughness interaction. 
The electron-electron interaction is taken into account by the Hartree potential in the single electron Hamiltonian. 
Details of the scattering mechanisms are presented in Appendix. 



VII. NUMERICAL MODEL 



The system under consideration consists of regions with high (the source and drain) and low (the channel) concentra- 
tions of electrons. The difference of the electron distribution in the corresponding regions would produce a considerable 
inaccuracy if we would have attempted to directly construct a finite-difference analog of Eq. ((19(1 . It is worth men- 



tioning that, in the quasi-classical limit, i.e. Ej S (( + rj/2) — Ej S (( — rj/2) 



Eq. 119(1 leads to the Boltzmann 



equation with an effective potential which has the following exact solution in equilibrium (Ep = Eps = Epp,)'- 



21.2 



exp [3 



K 2 k 
2mh 



+ E js (0 



2^,2 



hrp 

2m,yy 




(29) 
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For numerical calculations it is useful to write down the partial Wigner distribution function as fj S ((, k) — fj^(C> k) + 
ff s (C,k). Inserting this into Eq. H19(l . one obtains the following equation for / d s (£, fc) : 



W js (z,k,k')ff s (C,k')dk' = B js (C,k) , 



(30) 



where 



Bjs({,k) 



1 

2^ 



j +OC 

dk' J d v \E js (C 



dE JS (C ) 



V \sm[(k-k') V ]f^(Ck') 



(31) 



The unknown function fj S (C, k) takes values of the same order throughout the whole system, and therefore is suitable 
for numerical computations. In the present work, we have used the finite-difference model, which is described in 
Refi22iSi. The position variable takes the set of discrete values Q = A(i for {i = 0, . . . , Nr }. The values of k are also 
restricted to the discrete set k p = (2p — N k — 1) Ak/2 for {p = 1, . . . , N k }. The mesh spacing in the k space is 



Aj 



N k A( 



(32) 



The choice of discrete values for k follows from a desire to avoid the point k = and the need to satisfy a Fourier 
completeness relation. On a discrete mesh, the first derivative ^jp-(Ci) k p ) is approximated by the left-hand difference 
for k p > and the right-hand difference for k p < 0. It was shown in Ref^, that such a choice of the finite-difference 
representation for the derivatives leads to a stable discrete model. Projecting the equation H3()(l onto the finite- 
difference basis gives a matrix equation L • f = b. In the matrix L, only the diagonal blocks and one upper and one 
lower co-diagonal blocks are nonzero : 





/ 




-u 


. 





\ 




-V 


A 2 


-u . 





L = 







-V 


A 3 . 









V 








. 


■ A Nc 


-J 



(33) 



Here, A- U, and V are Nk x Nh matrices, 



2mfAC 



h 2 (2p-N k - l)Afc 



p" 



2mf AC 



h 2 (2p~ N k 



DAk ^ js ^ i: ^ p 



hP JS {k p ,k p/ )), (34) 



W\^=S^Bi^--j>Y lVU=Sp P >0[p-^-\, (35) 

and the vectors are 

[fi]p = fjs(Ci>kp), and [bi] p = B js (Q, k p ), i = 1, iVf_i, p = 1, N k . (36) 

A recursive algorithm is used to solve the matrix equation L • f = b. Invoking downward elimination, we are dealing 
with Bi = (Ai — V Bi-i)^ 1 !] and Cj — (Ai — VBi-i)~ 1 (bi + Vci^i) (i = 1, . . . , N^) as the relevant matrices and 
vectors. Then, upward elimination eventually yields the solution /; = Bif i+1 + a (i = — 1, . . . , 1). If a particular 
index of a matrix or a vector is smaller than 1 or larger than — 1, the corresponding term is supposed to vanish. 

In the channel, the difference between effective potentials Ej s (Q with different (j, s) is of the order of or larger than 
the thermal energy k^T. Therefore, in the channel only the lowest few inversion layer subbands must be taken into 
account. In the source and drain, however, many quantum states (J, s) of the motion along the x-axis are strongly 
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populated by electrons. Therefore, we should account for all of them in order to calculate the charge distribution. 
Here, we can use the fact that, according to our approximation, the current flows only through the lowest subbands 
in the channel. Hence, only for these subbands the partial Wigner function of electrons is non-equilibrium. In other 
subbands, electrons are maintained in the state of equilibrium, even when a bias is applied. So, in Eq. (|16|) for the 
electron density, we can substitute functions fj psps (z, k) of higher subbands by corresponding equilibrium functions. 
Formally, adding and subtracting the equilibrium functions for the lowest subbands in Eq. 1)16(1 . we arrive at the 
following equation for the electron density 

+00 

»(r) = n cq (r) + J- £ / dk (f js (z,k) |^ s (x;z)| 2 - f%(z, k) \4>g(x; z)f) , (37) 
y J' s -00 

where n cq (r) and ip^(x; z) are the electron density and the wavefunction of the radial motion in the state of equilibrium 
respectively. The summation in the right-hand side of Eq. is performed only over the lowest subbands. Since the 
electrostatic potential does not penetrate into the source and drain, we assume that the equilibrium electron density 
in these regions is described to sufficient accuracy by the Thomas-Fermi approximation : 

n cq (r) = N c ^=F 1/2 (/3(eV(r,z)+E F - E c )), (38) 

where the Fermi integral is defined by 

f°° ft 

exp( t -x) + l dt - < 3S » 

Here Ac is the effective density of states in the conduction band and Ep denotes the Fermi level of the system in 
equilibrium. 



VIII. RESULTS AND DISCUSSION 



For the numerical simulation of the SOI MOSFET we have fixed the geometrical parameters of the system as follows. 
The thicknesses of the semiconductor and oxide layers are chosen to be tsi = 20 nm and i ox = 2 nm respectively; the 
width of the MOSFET along the y-axis is taken to be 50 nm, the lengths of the drain and source are 20 nm each, 
while the barrier potential Vb at the Si/Si02 interface equals 3.2 eV. The values of the channel length considered 
for the numerical calculations are Lch = 20, 30, 40 nm. We further assume that the source electrode is grounded, 
whereas the potentials at the drain and gate electrodes equal Vb and Vg respectively. All numerical simulations are 
carried out with N$ = 240 and Nk = 120. The I - V characteristics are computed for a gate voltage Vq = 1.0 V, 
while the drain voltage Vb varies from to 1.0 V. The temperature of the system is T = 300 K. The doping level 
of the source and drain regions is Ad = 10 20 cm -3 , whereas the channel is doped with an acceptor concentration 
Aa = 10 18 cm -3 . In calculations of the scattering the following values of parameters for Si are used: A = 0.2, 0.4 nm, 
A = 1.5 nm 29 i 37 i 38 , E u =9.2 eV, p = 2.3283 • 10 3 kg/m 3 , v = 8.43 • 10 5 cm/ A 

The simulation of the SOI MOSFET requires a self-consistent solution of three equations: (i) the Poisson equation 
(0); (ii) the Schrodinger equation yielding the wavefunctions describing the transverse motion in the re-direction (0; 
(iii) the equation obeyed by the partial Wigner distribution function (|19|l . In the present calculations the eight lowest 
subbands are taken into account. 

Fig. 2 illustrates the distribution of the electrostatic potential for Vb = 0.5 V and for the value of the channel 
length Lqh = 30 nm. Clearly, the main part of the applied gate voltage voltage falls across the insulator while, 
due to the surface potential confinement, the current mainly flows in a thin layer near the Si/Si02 interface. At the 
p-n junctions (source-channel and drain-channel) the electrons meet barriers across the whole semiconductor. These 
barriers persist even for high values of the applied voltage Vb . 

In Fig. 3 the effective potential for the lowest inversion subband (j = 1, s = 1) is plotted as a function of z 
for different values of the applied bias ranging from to 0.5 V, Vg = 1 V, Lqh = 30 nm. In the case of ballistic 
transport (solid lines), the applied source-drain voltage sharply drops near the drain-channel junction. The scattering 
of electrons (dashed lines) smooths out the applied voltage, which is now varying almost linearly along the whole 
channel. The potential obtained by taking into account scattering is always higher than that of ballistic case. This 
can be explained by noticing that channel scattering locally increases the electron density and therefore enhances its 
capability of screening the applied gate voltage. 
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Fig. 4 shows the effective potential for the lowest inversion subband as a function of z for different values of the 
applied gate voltage Vq ranging from to 1.0 V, Vd = 0.5 V, L cii = 30 nm. The effective potential drop near the 
source-channel junction increases with increase gate voltage. This effect results into the potential barrier narrowing 
and lowering with the gate voltage. 

Fig. 5 shows the current density versus the source-drain voltage Vd for structures with channel lengths Lch = 20, 40 
nm. The solid lines represent the ballistic regime, the dashed and dotted lines reveal the effect of scattering due to 
ionized impurities, acoustic phonons and surface roughness at the Si/Si02 interface with an average displacement 
A = 0.2 nm and A = 0.4 nm respectively. At Vd = 0.2 V we observe a kink in the I - V characteristic. For Vb < 0.2 
V the derivative Vd/I reflects the differential resistance of the structure. Obviously, the inclusion of the electron 
scattering mechanisms enhances the resistance of the structure and softens the kink in the I - V characteristic, 
thereby appreciably degrading the current in the MOSFET. The collisions play an especially important role for low 
source-drain voltages. The surface roughness at the Si/SiC>2 interface gives rise to an appreciable though not dominant 
contribution to the scattering in the channel. When Vd exceeds 0.2 V, a "saturation regime" is attained although 
the currents do not converge towards constant values. From Fig. 5 one may conclude that a decrease of the channel 
length leads to an increase of the I - V slope in the "saturation regime" . This result is a consequence of a reduction 
of the p-n junction barrier potential as the length of the channel becomes shorter. 

Fig. 6 shows the I - V characteristics of the SOI MOSFET as a function of the gate voltage Vq. Fig. 7 illustrates 
the drain current I-q versus the gate voltage Vq of 30 nm channel-length SOI MOSFET for drain voltage Vd = 1 
V. We can see that the "threshold voltage" of the transistor is ~ 0.0 V. The influence of the scattering on the 
current diminishes as the the gate voltage lowers. This fact is due to a reduced electron concentration resulting into 
a lower amplitude of the scattering process. The decrease of the electron concentration, together with an increasing 
source-channel junction barrier, leads to the reduction of the current through the channel. Fig. 7 also informs us that 
we are still far removed from the International Technology Roadmap for Semiconductors (ITRS) requirements 40 . In 
particular, whereas the on-state current (~ 5 A/cm) is in the range of the requirements (~ 7.5 A/cm), the off-state 
current (0.5 A/cm) is far above target (~ 0.0015 A/cm). Therefore, further optimization of the silicon-film thickness 
as well as channel-doping engineering is needed. 

In Fig. 8 the contours of the absolute value of the ground-state Wigner function (/n) are plotted with and without 
the presence of scattering for structures with channel lengths Lqh — 20 nm (the first row), 30 nm (the second 
row), 40 nm (the third row). In these figures the first column corresponds to the ballistic regime, the second and 
the third correspond to the regime with scattering due to ionized impurities, acoustic phonons and surface roughness 
at the Si/Si02 interface with the average displacement A = 0.2 nm and A = 0.4 nm respectively. Darker areas 
in these plots indicate higher density of electrons. Far away from the p-n junctions where the effective potential is 
varying almost linearly with z, the partial Wigner function is positive definite and can be interpreted as a probability 
distribution of electrons in phase space. When electrons are accelerated in the inversion layer without scattering, 
their velocity increases monotonously along the whole channel. Therefore, in the phase-space representation, the 
distribution of ballistic electrons looks as a narrow jet stream pointing from the source to the drain along the channel, 
while scattering is seen to wash out the electron jet in the channel. Consequently, the electron transport through the 
channel combines the elements of both diffusive and ballistic motion. The mean-free-path less than the device length 
does not mean that the device is near-ballistic. 

In this paper we have studied in detail the effects of ultra-short channels in an SOI MOSFET. To benefit from the 
SOI architecture, the thickness of the channel should be reduced to 10 nm or less. The simultaneous scaling of the 
channel length and thickness is currently under study. 
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APPENDIX: SCATTERING MECHANISMS 

We consider the electron scattering due to impurities, acoustic phonons and surface roughness at the Si-Si02 
interface. In case of the electron scattering due to impurities 34,41 H ln t is the electron-impurity potential minus its 
spatial average. The latter term has already been included in energy level calculations (band bending calculations), 
when solving Poisson and Wigner equations self-consistently. This spatial average being, by definition, translation 
invariant in the layer plane, has a zero matrix element between initial and final plane wave states and can thus be 
dropped. The impurities are assumed to be randomly placed at the sites Hi = (pi,Zi). The Hamiltonian describing 
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scattering due to the impurities in the inversion layer is given by 

flint 0) =J2Vq (zj exp (iq (p - pj) . 



(A.l) 



where 



2^e 2 



ei (g + g s F(g)) 



(A.2) 



Here g s is the Thomas-Fermi screening constant and F (q) and _F (g, Zi) are form-factors, corresponding to the screening 
impurity potential^. Because the screening effect in the channel is very large, as seen in the Refsi22*2i, we can model 
the scattering potential of an ionized impurity as U(r) = (Aire 2 R 2 / 'e i)<5 (r) , where R s determines a cross-section for 
scattering by an impurity. Consequently, the absolute value of the matrix element is 



\{jsp'k'\ U(r - Ri) \js P k)\ = A-Ke'Ri/e 1 ^ js {x i ,z i ). 
Averaging this over a uniform distribution of impurities results in the following scattering rate 

p ™k,] S p>ki 0) = 27tCim / dxipj s (x, z)8 (E jB p> k > - E jspk ) , 



(A.3) 



(A.4) 



where 



ClM 



1 „ /4ie 2 iif 2 



hS V £1 



(A.5) 



N\ is the acceptor concentration and S = L y L z . 

Lattice vibrations are an inevitable source of scattering, and we need to take them into account when we calculate 
device characteristics at room temperature. In this work we restrict ourselves to intra-valley and intra-subband 
acoustic phonon scattering. In the approximation of an isotropic elastic continuum we have 



rVe-ph _ ~ 
•"int — 



d 

D V • u(p,z) + -g-u z (p, z) 



q 



(7 q 6 q exp(iq • r) + h.c.) , 



(A.6) 



within 



bq 



pVLOq 



EE, p and o; q respectively denote the deformation potential, the mass density of the semiconductor and phonon 
frequency. The matrix elements describing single phonon emission and absorption processes are respectively given by 



jsp'k';Na + l 



fl 



e— ph 



jspk]Nq) = TqV^q + 1 U s P' k ' |exp(-iq • r)\jspk) 



(A.7) 



jsp'k'-N^-l 



fl. 



c— ph 



jspk;Nq) = ~fq^/Nq(jsp'k' |exp(iq • v)\jspk) 



where A q is a phonon occupation number. The emission (+) and absorption (-) rates are 



P 



jspk,jsp' k' 



2tt 



p jspk,js P > w = ~r l 7 ql 2 N i I l ex P( ic t 1 r ) 1 3spk)\ 2 S {E Jsp , k , ~ E jspk - hu q 



X! lTq| 2 (A q + 1) \ (jsp'k' |exp(-iq • r)| jspk)f 8 (E jsp r k > - E jspk + fuv q 



(A.8) 

(A.9) 
(A.10) 



At room temperature both expressions are approximately the same. Indeed, for T — 300 K the thermal energy k^T 
largely exceeds the phonon energy fiw q , and therefore the acoustic deformation can be considered elastic. For low 
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energies, we can approximate the phonon occupation number as Nq ~ k^T/hujq 3> 1 and adopt the Debye model 
w q = v oQ where Vq is the sound velocity. Assuming equipartition of energy in the acoustic modes, we obtain the 
scattering rate : 



Pj'apkjsp'k' 0) = SttCph / dxrfjx, z) 5 (E jBp > k > - E jspk ) . 



(A.ll) 



with 



PH 



1 „ 2 k B T 
HS^ pv 2 ' 



(A.12) 



Ideally, the Si-Si02 interfaces would coincide with a planes x — and x — igi- The real interfaces however are known 
to fluctuate randomly around these planes and, under the assumption that the boundary between Si and Si-Si02 is 
still abrupt, the local deviations due to surface roughness may be characterized by a quasi-continuous function A(p). 
Here p = (y, z) represents the two-dimensional position vector parallel with the ideal interfaces. We further assume 
that the surface potential varies linearly with A, i.e. 



V (x + A (p)) pa V (x) + A (p) 



dV jx) 
dx 



(A.13) 



Neglecting the effect of surface roughness on the subband energies and wave functions, we may calculate the scattering 
matrix element for surface roughness to the lowest order in A as follows 29 - 43 : 



j S pk\H?*\jsp'k' 



dipjs dip. 



2m* x dx dx 



*k-k' 



*k-k' 



a:=0 



da; | ipj 



dV (x) 
dx 



where k = and Ak— k' is the Fourier transform of the A (p) 

A(p) = ^A q exp(iq- p) 



(A.14) 



(A.15) 



In the Born approximation, only the magnitude squared of Ak_k' (referred to as the power spectrum) is needed, and 
thus the phase of Ak_k' can be neglected. We propose a Gaussian distribution to describe the spatial correlation of 
the surface roughness 



C(p-p') = (A(p)A(p'))=A 2 exp - 



\P-P 

A 2 



(A.16) 



where A is the root-mean-square value of A(p), and A, is referred to as the correlation length governing the the decay 
of the auto-correlation function. In the above expression the brackets ( ) denote sample averages 29 . By convolution, 
the power spectrum reduces to the Fourier transform of the auto-corrleation function l|A.16(l , that is given by 



5(q) - |A q | 2 = 7rA 2 A 2 exp (-q 2 A 2 /4) 
and where q is the scattered wave vector. 

Consequently, the interface roughness scattering rate takes the form : 



(A.17) 



pSR ( \ 

jspk,jsp' k' \ J 



2tt 2 



A 2 A 2 exp 



hS 

X $ (Ejsm'k 



|k-k'l 2 A 2 



dx\ipj S (x,z)\' 



E H 



The total scattering rate 



P , ■ ,u - P IM 
r jspk,]sp'k' — r j S pk,jsp'k' 



pPH 

jspk.jsp' k' 



pSR 

j spk,j sp' k' 



is then inserted into Eq. (|26|l in order to obtain the one-dimensional scattering rate 



dV(x) 
dx 



(A.18) 
(A.19) 



P 



IM+PH 



(z, k,k ) = m 



2L 2 V (3 



h 2 TT 



(Cim + C P h)F /3 



h 2 k 2 
2m? 



h 2 k r 
2m? 



tpj s (x,z) dx, 



(A.20) 
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2a2 ^yy y 



i^ R (z,M') = ^A 2 A 2 (m 



1 H 2 TT 



where 



exp 



(k - k') 2 A 2 



n (k 2 - k r 



Ax\ip js {x,z)\ 



dV (x) 
dx 



(A.21) 



n(aO 



exp 



dp- 



h 2 p 2 



2y^ 



Vp 2 + ■ 



( 



exp 



(p - y/p 2 + xj A 2 



exp 



(p + \Jp 2 + x\ A 2 



X 




(A.22) 



and K (x) is a McDonald function 44 . 



F(x) =- x ' 2 X (|z|/2), (A.23) 
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FIG. 1: Scheme of the SOI MOSFET. 




FIG. 2: The electrostatic potential energy in the SOI MOSFET with Lew. = 40 nm at V G = 1 V and Vu = 0.5 V. 




FIG. 3: Effective potential as a function of z for various Vb, £ch = 30 nm. 




FIG. 4: Effective potential as a function of z for various Vg, £ch = 30 nm. 




FIG. 5: Current- voltage characteristics at Vg = 1 V for different channel lengths. 
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FIG. 6: Current-voltage characteristics for SOI MOSFET with Lew. = 30 nm. 




FIG. 7: Drain current vs. gate voltage characteristics for drain voltage Vd = 1 V. 
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FIG. 8: Contour plots of the absolute value of the partitial Wigner distribution function |/n (z, k)\ at Vg = 1 V, Vb = 0.5 V 
and various channel length. Leu = 20 nm: (a) without scattering, (b) with scattering A = 0.2 nm, (c) with scattering A = 0.4 
nm; Leu = 30 nm: (d) without scattering, (e) with scattering A = 0.2 nm, (f) with scattering A = 0.4 nm; Leu = 40 nm: (g) 
without scattering, (h) with scattering A = 0.2 nm, (i) with scattering A = 0.4 nm. 



